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GLOBAL  RECEIVER  OPERATING  CHARACTERISTIC  (ROC)  CURVES 
FOR  THREE  NETWORK  TOPOLOGIES 

1.  INTRODUCTION 


In  traditional  multisensor  detection,  local  sensors  transmit  data  to  a  central  processor  that 
performs  optimal  detection  based  on  conventional  statistical  techniques.  In  decentralized  or 
distributed  detection,  each  geographically  dispersed  sensor  executes  limited  data  processing 
before  transmitting  condensed  information  to  a  central  site.  Thus,  the  sensor  becomes  a 
processor.  Interest  in  distributed  detection  has  intensified  due  to  the  relatively  low  cost  of 
sensors,  the  inherent  network  redundancy,  the  availability  of  high-speed  communication 
networks,  and  the  increased  computational  capacity  of  sensors.1,2 

One  element  of  distributed  detection  is  receiver  operating  characteristic  (ROC)  curves, 
which  show  the  relationship  between  the  probability  of  detection  and  the  probability  of  a  false 
alarm.  This  document  surveys  the  results  of  multisensor  detection  work,  with  a  focus  on  how  to 
combine  multiple  ROC  curves  into  a  single  network  or  global  ROC  curve  to  describe  the 
network.  The  solution  depends  on  the  network  topology.  While  other  topologies  exist,  this 
document  examines  only  three  networks: 

•  a  fusion  network  (e.g.,  common  undersea  picture), 

•  a  serial  network  (e.g.,  large  numbers  of  distributed  sensors),  and 

•  a  tree  network  (e.g.,  command  and  control). 

The  fusion  and  serial  networks  are  special  cases  of  the  tree  network.  For  each  of  the  three 
topologies  considered,  formulas  are  presented  to  compute  the  global  ROC  curve.  The  focus  is  on 
the  fusion  network,  however,  with  minimal  results  presented  for  the  serial  and  tree  networks. 

For  each  topology,  the  network  has  N  local  processors  conducting  binary  hypothesis  testing. 
Processor  i  receives  input,  decides  between  two  hypotheses  (Ho  and  H\),  and  outputs  a  decision 
Uj,  which  is  determined  as  follows: 

f  0  processor  i  decides  H0 

ut  =  \  .  (1) 

[1  processor  /  decides  Hx 

The  processors  function  together  to  generate  the  best  network  decision  D.  Thus,  the  network 
decision  follows  an  optimality  criterion.  The  Neyman-Pearson  optimality  criterion,  used  in  this 
analysis,  constrains  the  global  probability  of  false  alarm  PF<  S  and  maximizes  the  global 
probability  of  detection  Pd  using  Lagrangian  multipliers.  The  Lagrangian  function  is 
P(PF  -S)-PD,  where  p  is  the  Lagrangian  multiplier  or  the  global  threshold  that  yields  the 

maximum  network  probability  of  false  alarm  such  that  PF<  8.  The  Bayesian  minimum  cost 
criterion  is  addressed  in  section  5,  “Bayesian  Formulation.” 
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For  this  discussion,  the  decision-making  processes  of  the  local  processors  are  assumed  to  be 
independent.  Associated  with  each  of  the  N  local  processors  is  a  ROC  curve.  For  each  processor 
i,  the  probability  of  false  alarm  and  the  probability  of  detection  (independent  of  network 
topology)  are  defined  as3 


oo 


Pf  =  jP(x\H0)dx  =  Pr(Ui=\\H0), 

ai 

(2a) 

oo 

Pdi  =  \P(x \H,)dx -Priu.-UH,), 

(2b) 

where  a,  is  the  detection  threshold,  Pr(itj  =  1  |  Ho)  is  the  probability  of  deciding  H\  is  true  when 
Hq  is  true,  and  Pr(u,  =  1  |  H\)  is  the  probability  of  deciding  H\  is  true  when  H\  is  true. 
(Throughout  this  document,  a  lowercase  subscript  refers  to  an  individual  sensor,  while  an 
uppercase  subscript  refers  to  the  network.)  The  definition  of  the  detection  threshold  will  depend 
on  the  problem  being  analyzed.  For  example,  for  a  sonar  detection-related  problem,  the  detection 
threshold  might  correspond  to  the  signal-to-noise  ratio  required  for  detection. 

An  issue  arises  at  the  fusion  center  when  multiple  processors  decide  on  the  same  hypothesis 
when  the  problem  assumptions  (known  at  the  fusion  center)  prevent  this  hypothesis  from  being 
correct.  For  example,  a  single  red  submarine  is  in  a  known  area  (area  search),  and  N  individual 
processors  are  tasked  to  detect  the  submarine  in  non-overlapping  sub-regions  over  some  period. 
How  the  individual  decisions  are  combined  into  one  solution  resolves  the  issue  of  the  network 
decision  when  multiple  processors  detect  a  submarine  at  the  same  time. 
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2.  FUSION  NETWORK 


2.1  FUSION  CENTER  DOES  NOT  RECEIVE  EXTERNAL  OBSERVATION 

Figure  1  shows  the  topology  of  the  JV-processor  fusion  network  where  the  fusion  center 
receives  only  the  output  of  the  individual  local  processors.  Each  local  processor  receives  an 
external  observation  vector  Xh  and  outputs  a  decision  The  fusion  center  receives  the  N 
decisions  ui,...,Un,  combines  the  individual  decisions  according  to  some  fusion  rule,  and  outputs 
a  final  network  decision  D. 


Figure  1.  N-Processor  Fusion  Network 


The  general  fusion  rule  is  given  by4 
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Pi*,..  ■?%)  —  ] 


0 


if  1  ^l)  V,  p 

P{ux,u2 uN  |  H0 ) 
if  1  H\)  <  p 

P(liX  ,^2,  .  .  .  5  Ufl  |Ffg) 


(3) 


where  f3  is  the  global  threshold,  which  is  the  Lagrangian  multiplier  for  the  Neyman-Pearson 
criterion.4  Since  the  local  decisions  are  independent,  the  likelihood  ratio  can  be  simplified  to4 
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.  Ht  , . . . ,  |  //l:  )  i  P(»,  |  H ,  I 

nj’(",=l|g1)rff(-',=0|W,) 

i<ip(li,=l|ffo)i„1)P(»,=0|ff1>) 


=nfn 

t/(l)  V,  C/(0) 


1-p 


4 


1-p, 


(4) 


where  (7(1)  is  the  set  of  all  i  such  that  w,  =  1,  (7(0)  is  the  set  of  all  i  such  that  m,  =  0,  and  Pd  and 
Pf  represent  the  probability  of  detection  and  the  probability  of  false  alarm  for  processor  /, 
respectively. 

The  (Pf,  Pd)  operating  point  of  the  fusion  center  is4 
PF  —  ^ F(u j , . . . , uN )P(u1 uN  \  H0 ) 

«1  “2  «N 

ax  u2  uN  i=l 

. <! 

«i  u2  uN  U(  1)  C/(0) 


and 


N^P(.U\im  •  I  ^l) 

«1  u2  uN 

'ZZ-z^.-.^tn^if/,) 

«1  «2  1=1 

■zz-z^.-.«.)n^no-^)-  (5b> 

«,  «2  C/(l)  L?(0) 

Equations  (5a)  and  (5b)  can  be  difficult  to  compute  for  large  N. 

There  is  one  caveat  when  using  Lagrangian  multipliers:  the  individual  ROC  curves  must  be 
convex.  Otherwise,  the  solution  might  not  be  locally  optimal. 

A  special  case  that  can  be  efficiently  evaluated  is  when  the  fusion  rule  is  a  Ai-of-out-./V  logical 
decision  (i.e.,  k  processors  decide  hypothesis  H\).  If  the  operating  points  of  the  local  processors 
are  identical,  then  given  a  point  (Pf,  PJ)  on  the  individual  ROC  curve,  equations  (5a)  and  (5b) 
determine  a  point  on  the  global  ROC  curve  as4 
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(6a) 


k- 1 


pF= i-X 

j= 0 


"TV" 

J  > 


<p,y(i -/»,)" 


and 


i-1 


^=i-Z 

7= 0 


'N' 
J  y 


\(pdy(i-pd)N  j.  (6b) 

The  global  performance  can  also  be  expressed  in  terms  of  the  incomplete  beta  function5 

Ix(a,b)  =  — - —  \ta\\-t)b-ldt 
B(a,b)  * 

X 

j>  \l-t)bldt 


_  0 


\ta\\-t)bldt 

0 

where  B(a,b )  is  the  complete  beta  function.  Thus,  the  global  performance  is4 
pF=iPa,N-k+ 1) 


(7) 


and 


PD=IP(k,N-k+ 1). 


(8a) 


(8b) 


If  Py<  0.5  (which  is  almost  certainly  true  in  applications),  then  N  >  2.  Using  equations  (8a) 
and  (8b),  the  global  probability  of  false  alarm  and  global  probability  of  detection  can  be 
expressed  in  terms  of  the  approximation  shown  in  appendix  A  and  the  optimal  fusion  rule  A:4’5 


PF{k)  =  Q\ 


and 


PD(k)  =  Q\ 


(l -  Pf  f*  (9k  -  l\N  -  k  +  if3  -  Pf%  {9(N  -  k  + 1) - 
3  k^jN-  k  +  l^p/^k'^ +(l-PfY(N-k  +  l)^ 


(l -PdY (9 k - 1 \N-k  +  lft  - p/> (9 (TV -k  + 1)- 1  )k% 
3k%  y/N-k  +  l^p/ik^ +(l-PdY(N-k  +  lY  j 


(9a) 


(9b) 
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Equations  (9a)  and  (9b)  will  produce  the  optimal  solution  if  k  and  a  are  the  optimal  solutions  to 
the  global  threshold  equations. 


2.2  FAULT-TOLERANT  REDUNDANT  NETWORK 

Figure  2  shows  a  fault-tolerant  redundant  fusion  network.  Each  of  the  N  identical  local 
processors  receives  the  same  observation  vector  X,  and  outputs  a  decision  If  no  faults  occur, 
the  processors  output  identical  decisions.  The  fusion  center  combines  the  N  decisions  uj,  ...  uN  to 
reduce  the  effect  of  a  fault  in  any  local  processor,  and  the  fusion  center  outputs  a  final  decision 
D. 


Figure  2.  N-Processor  Redundant  Fusion  Network 


Network  faults  can  be  classified  into  one  of  two  types:  hard  failures,  in  which  the  network, 
or  one  of  its  elements,  is  not  able  to  deliver  any  traffic  at  all;  and  soft  failures,  which  includes 
network/service  anomalies  or  degradations  in  various  performance  parameters.  (The  cause  of  the 
fault  is  not  significant,  and  therefore  is  not  discussed  here.)  A  failure  occurs  when  the  processor 
does  not  deliver  the  theoretically  expected  decision.  Since  the  fusion  processor  is  simpler  than 
the  local  processors,  the  probability  of  failure  in  the  fusion  processor  is  insignificant  compared  to 
the  probability  of  failure  in  the  local  processor.  Therefore,  the  only  faults  considered  in  this 
paper  are  those  located  at  the  individual  processors,  and  not  at  the  fusion  center. 
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Processor  performance  is  expressed  as  a  function  of  the  fault-free  performance  (i.e.,  the 
theoretically  optimal  value).  Assume  that  these  probabilities  are  independent  of  the  hypothesis 
and  observation.  For  processor  i,  the  probability  of  false  alarm  is4,6 

Pf=P(u,=\\Ha) 

=  P{ut  =  1 1 M,  =  l)p(s;  =l|//0)+  p(ut  =  1 1  ut=  0)p(k,  =  0|  tf0) 

=  (l-Po)Pfl+Pi^-Pfl) 

=  Pi+(l~Pi-Po)Pf>  (] 


where 

Uj  =  the  observed  decision, 

ui  =  the  theoretically  optimum  decision, 

the  superscript  ~  indicates  that  the  value  is  the  theoretically  optimal  value, 

P(u:  =  1  j  u*  =  0)  =  p1 ,  and 

P(M.  =  0 1  u*  =  1)  =  p0 . 

Similarly,  the  probability  of  detection  is4,6 

pd=pM^-p,-Po)pdr  (ii) 

Note  that  the  original  ROC  curve  for  processor  i  is  scaled  by  a  factor  1  -  p\  -  po  in  each 
dimension. 

When  the  fusion  rule  is  a  £-out-of-jV  logical  function,  the  probability  of  false  alarm  for  the 
redundant  network  can  be  written  as4 

PP=P{D  =  1\H0) 


( 

' N ^ 

\ 

/  ~  \ »  /  _ _ _  vy  j 

V>  =  l|tfo) 

1 

M 

P{u,  =1\D  =  l,//0 //>(»,  =  0 1 D  =  l,Hj 

^  j= 0  1 

J  J 

/ 

j 

(  u-i 

/  »v  \i  (  ~  XN  i  1 

o 

II 

1-1 

p(u,  =  1 1 D  =  O'Hjpfa,  =  0 1 D  =  0,HJ 

i 

l  J* 

J  ) 

J 

(  »-i 

(N') 

> 

1  /  vl 

f  k- 1 

'N' 

\ 

PF 

1-Z 

(i  -PoyPo~j 

< 

1 

+ 

i-Z 

(i -pjpr1 

l  j=0 

J  j 

j 

i 

l  >• 

J  J 

J 

Similarly,  the  probability  of  detection  can  be  written  as4 


(12) 
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2.3  FUSION  NETWORK  WITH  FEEDBACK 

Figure  3  shows  a  fusion  network  with  feedback.  On  every  time  step,  the  N  local  processors 
communicate  their  decisions  to  the  fusion  center,  and  the  fusion  center  then  communicates  the 
global  decision  back  to  the  N  local  processors.  The  system  operates  as  follows: 

1.  At  time  step  t,  the  k .  processor  makes  a  local  decision  u[ ,  based  on  the  previous  global 
decision  D‘  \  the  current  observation,  and  all  previous  observations. 

2.  The  fusion  center  receives  all  local  decisions,  and  combines  them  to  generate  the  global 
decision  Df. 

3.  The  global  decision  is  then  transmitted  to  all  local  processors  for  use  at  time  step  t  +  1 . 


Figure  3.  N-Processor  Feedback  Network 

The  A  local  processors  generate  a  set  of  decisions  At  =  {u[,ut2,...,u‘N),  t=  1,2,...  in  the  tth 
processing  interval.  The  binary  decisions  in  At  regard  the  hypotheses  Ho  and  H\.  The  set  A,  is 
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combined  at  the  fusion  center  to  form  D'.  From  the  discrete  version  of  the  Neyman-Pearson 

•  •  ■  7  ^ 

lemma,  the  most  powerful  combination  rule  at  the  fusion  center  is 


4>,U)= 


P(A,  |  //,)>  XIP(AI  I  Ha) 


< 


(14) 


where  y  e  [0,1],  and  Xt  >  0  is  the  Lagrange  multiplier. 

At  time  step  t,  the  probability  of  false  alarm  can  be  written  as7 


P,F^E{ct>l{A,)\Hy}  = 


A 

+ 0  -  K  '  )£  *<  (A  )p(A- 1  D"‘  =  °.  ■ ", ) 

A 

YMA)p(A\h„) 

.  A 

•  •  •  *7 

Similarly,  the  probability  of  detection  can  be  written  as 


t>  2 

t  =  l 


(15) 


4 


P'=2?M(4)|tf0}  = 


+(i-C')Z«s.(4)/'(4I£>m=o,//,) 

4 


t  >  2  . 


t  =  l 


(16) 


An  alternate  topology  for  a  fusion  network  with  feedback  uses  parleying.  This  organization 
operates  as  follows:8 

1 .  Each  sensor  receives  an  observation  vector  Xj. 

2.  Each  sensor  makes  a  tentative  decision  based  on  its  own  data,  and  transmits  this  decision 
to  all  the  other  sensors  in  the  network. 

3.  Based  on  its  own  original  data  and  the  tentative  decisions  of  all  the  other  sensors,  each 
sensor  then  “rethinks”  and  possibly  revises  its  previous  decision.  This  decision  is  then  broadcast 
to  all  the  other  sensors  in  the  network. 

4.  The  parleying  in  step  3  continues  until  all  sensors  agree. 

Although  this  topology  has  application  in  military  analysis,  it  is  not  addressed  in  this  document. 
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3.  SERIAL  NETWORK 


Figure  4  shows  an  TV-processor  serial  network.  The  first  processor  receives  only  an  external 
observation  vector  X\  and  then  outputs  its  decision  u\.  The  remaining  processors  receive  an 
external  observation  vector  Xt  and  a  decision  Uj.\,  and  then  according  to  some  fusion  rule  output 
decision  w,.  Processor  N  outputs  the  decision  of  the  network.  There  are  two  main  reliability 
problems  with  a  serial  network.2  The  first,  and  most  serious  issue  of  the  two,  is  communication 
failure.  If  a  communication  link  fails  anywhere  in  the  chain,  then  all  the  processors  before  the 
failure  are  effectively  removed  from  the  system,  resulting  in  serious  performance  degradation. 
The  second  problem  is  one  of  accumulated  delays  because  each  processor  has  to  wait  for  input 
from  the  previous  processors  before  generating  a  decision. 


X, 


X. 


Figure  4.  N-Processor  Serial  Network 


For  N  serial  processors,  the  probability  of  false  alarm  and  probability  of  detection  are 
computed  recursively.  Consider  a  serial  network  with  two  local  processors.  Processor  one 
outputs  decision  u\,  and  processor  two  outputs  decision  D.  The  network  probability  of  false 
alarm  is  determined  using  the  theorem  of  total  probability4 

Pf=P(D  =  1\H0) 

=  P(D  =  1 1  Mj  =  1, H0)P(Ul  =  1 1  H0)  +  P(D  =  1 1  u,  =  0, H0)P{Ul  =  0 1  H0).  (17) 

Replacing  Ho  with  H\  in  equation  (17)  gives  the  probability  of  detection.  The  optimal  thresholds 
for  each  processor  are  evaluated  using  the  probabilities  of  false  alarm  and  detection  for  the  local 
processors.  The  thresholds  for  each  processor  are  defined  as  aj ,  the  threshold  used  by  processor 

i  when  it  receives  decision  u,.\  =  j.  For  the  two-processor  serial  network,  the  complete  set  of 
threshold  equations  is  4 
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PfM)-Pf2(a02 

P  Pd2{a\)-Pd2(al) 


<A=P 


]zW> 

1  -PdlW 


<*\=P 


PflM 

PdM) 


Since  P(u2  =  1 1  ux  =  i,  H0 )  =  Pf  (a‘2 ) ,  the  global  probability  of  false  alarm  is4 


PF  =  p,t  (arj )PA(al)  +  ph  (a2°)(l -  J»  (a, )) . 
Similarly,  the  global  probability  of  detection  is4 

Pd  =  PdS^)Pd^)^Pd2{ali\-Pdi{aS). 


(18) 

(19a) 

(19b) 


Now  consider  an  //-processor  serial  network.  The  serial  network  with  j  processors  can  be 
considered  as  a  serial  network  of  two  processors.  The  first  “processor”  has  the  performance  of 
the  first  j  -  1  processors,  while  the  second  “processor”  is  the  7th  processor  of  the  network.  The 
performance  of  local  processor  i  is  Pd  and  Pf ,  while  the  performance  of  a  network  with  j 

processors  is  PDj  and  PF .  Recursive  equation  (20)  provides  the  solution  to  the  combined  ROC 
curve  for  the  network:4 


aj-i=Pj 


Pfj(a))-Pfj(a^) 

PdWj)-PAa]) 


?f  ( ai  1) 

Oo  »  fj-ik 

3  Pd^Mj-i) 

'  Pi\  -PD(aMy 


a 


pFj = p,„  0 a («;> + (1  -  p,„  >K  <«;>. 

PD,  =  Pd„  &J-M,  («;  )  +  (l  -  Pd,.,  )K  K  ) ' 


(20) 


for  j  =  2,. .  .JV,  where  PF  =  Pf> ,  PD]  =  Pd  ,  and  J3j  is  the  global  threshold  for  a  system  of j 

processors  computed  as  above.  Note  that  these  equations  are  coupled  and,  in  general,  are  quite 
difficult  to  solve. 
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If  the  ROC  curves  for  the  local  processors  are  symmetric,  then  the  threshold  equations  are 
greatly  simplified.4  Symmetric  local  ROC  curves  imply  that  the  global  ROC  curve  is  also 
symmetric.  Thus,  the  global  threshold  J3=  fa  =  1  implies  that  a,  =  1  for  j  =  l,2,...,yV-l .  In 
addition,  PD  =  1  -PF  for  j=  1,2,...^V.  As  a  result,  the  equations  in  (20)  become4 


ocM=l, 


1-P„ 


a°  = 


rj- 1 


rj- 1 


(21) 


for  j  =  2,.  ..JV,  where  PF  =  Pf  and  PDi  =  Pd  .  Equation  (21)  is  then  relatively  easy  to  solve. 
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4.  TREE  NETWORK 


A  tree  topology  is  more  difficult  to  analyze  than  either  the  fusion  or  the  serial  network. 
Figure  5  shows  a  distributed  detection  network  with  a  tree  structure.  Node  N  in  figure  5  outputs 
D,  the  final  decision  of  the  network.  A  tree  network  can  be  represented  by  a  directed  acyclic 
graph  with  each  node  being  a  local  processor.  Edges  represent  one-way  communication  links 
between  processors. 


Figure  5.  Distributed  Detection  Network  in  a  Tree  Structure 
(node  N  outputs  final  decision) 


Each  of  the  leaf  nodes  in  a  tree  structure  receives  an  observation  vector  A,,  makes  a  decision 
Uj  based  on  that  event,  and  transmits  its  decision  to  each  of  its  immediate  successor  nodes.  The 
internal  root  node  receives  decisions  from  all  its  predecessor  nodes,  combines  the  decisions 
according  to  some  fusion  rule,  and  transmits  the  decision  to  all  its  immediate  successor  nodes. 
For  example,  in  figure  5,  node  j  receives  observation  vector  Xj,  makes  a  decision  Uj,  and  transfers 
this  decision  to  node  /.  Node  i  receives  four  incoming  decisions,  combines  them  according  to 
some  fusion  rule,  and  transmits  decision  w,  to  its  successor  node. 

Pifi)  =  P(uj  =  1  \H0)  and  Po(i)  =  P(«,  =  1  \H\)  are,  respectively,  the  probability  of  false  alarm 
and  the  probability  of  detection  of  a  subnetwork  with  node  i  as  the  root  node.  Using  the  theorem 
of  total  probability,  if  D  is  the  final  decision  of  the  network,  the  probability  of  detection  and  the 
probability  of  false  alarm  can  be  written  (independent  of  the  optimality  criterion)  as9,10 
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(22a) 


PD  =  P{D  =  \\uk=0,H1) 

+  [P{D  =  \\uk=\,Hl)-P{D  =  \\uk=Q,Hl)]PD(ky 
and 


Pp=P(D  =  \\uk=Q,H„) 

+  [P(D  =  1 1  u„  =  P(D  =  1 1 «,  =  0,  tf0)]/V(*)’ 

for  any  node  k  where  k*N. 


(22b) 
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5.  BAYESIAN  FORMULATION 


Until  this  point,  the  network  was  assumed  to  follow  a  Neyman-Pearson  optimality  criterion. 
Now  assume  that  the  network  follows  a  Bayesian  minimum  cost  optimality  criterion.  The 
objective  of  the  Bayesian  criterion  is  to  minimize  the  overall  cost  of  the  network  decision.  If 
J(D,  Hj )  is  the  cost  of  the  final  processor  deciding  D  when  hypothesis  Hj  is  true,  optimality 
demands  that  the  total  expected  cost  E[J(D,  HJ)]  be  minimized.  The  functions  that  characterize 
the  local  processor  decisions  are  chosen  such  that  E[J(D,  HJ)]  is  a  global  minimum.  This 
formulation  could  be  difficult  to  implement  in  practice,  since  a  priori  the  probabilities  of  each 
hypothesis  and  the  costs  associated  with  false  alarms  and  detection  are  needed. 

The  only  changes  are  for  the  network  where  the  fusion  center  does  not  receive  an 
observation  vector,  and  the  fusion  center  network  with  feedback.  All  the  equations  for  the  other 
networks  remain  the  same.  For  the  fusion  center  network  where  the  fusion  center  does  not 
receive  an  external  observation  vector  (see  section  2.1),  the  only  change  is  to  /?,n  or 

P(ff,)[./<0, //,)-■/(!,«,)]  ' 

In  section  2.3,  the  global  probability  of  false  alarm  and  the  global  probability  of  detection 
change.  The  system  probability  of  error  is12 

7>;=p;p(//0)+p^(/f1>.  (24) 

Expanding  P'F  in  terms  of  Dt]  12 

P'F  =  P(D'  =  1 1  H0 )+  P(D'  =  1 1  D'1  =  1, H0 )p(d'  1  =  1 1 H0) 

+  p(Dt=l\Dl^O,H0)p(D,-l=0\H0)  .  (25) 

Replacing  p(d'1  =  0 1  H0 )  with  1  -  p(d'1  =  1 1  H0)  and  rearranging  terms  yields12 

P‘  =  p{d'~1  =  1 1  H0  \P(D‘  =  1 1 D'1  =  1,  H0 )-  P(D'  =  1 1  £>m  =  0,  H0 )] 

+  p(d'  =l\D‘  1  =0,Ho)  ■  (26) 

Thus,  the  global  probability  of  false  alarm  on  time  step  t  can  be  written  as12 

=  /■;-|[/’;(d-1  =  i)-  p'f(d'-' = o)]+  p’(d'-' = o)  •  (27) 

Similarly,  the  global  probability  of  detection  on  time  step  t  is12 
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/>k=oi//„) 

1  -Pfh(u'0'  =0 )-P'(u',-'  =l)]+7S(«'-‘  =l)  • 


6.  SUMMARY 


This  document  surveyed  work  on  how  to  combine  multiple  ROC  curves  into  a  single 
network  ROC  curve  for  three  network  topologies.  In  warfare  analysis,  obtaining  a  single 
descriptor  of  a  system  is  critical.  As  the  focus  on  networking  and  information  technology  in 
warfare  analysis  intensifies,  the  ability  to  describe  a  system  quickly  and  easily  will  be 
imperative.  A  network  ROC  curve  is  a  rapid  and  simple  way  to  explain  the  performance  of  a 
system.  While  many  of  the  formulas  presented  here  are  difficult  to  solve,  through  parametric 
analysis  they  can  provide  insight  into  trends  in  various  warfare  areas  where  ROC  curves  arise, 
including  detection/classification,  command  and  control,  and  decision-making. 
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APPENDIX  A 

APPROXIMATION  OF  INCOMPLETE  BETA  FUNCTION 


When  (a  +  b  -  1)(1  -  x)  >  0.8,  the  incomplete  beta  function  can  be  approximated  by  the 
complement  of  the  cumulative  Gaussian  distribution  function.5  If  a  +  b  >  6,  then5 


^  co  -t2 

Ix(a’b)  =  Q(y )  +  £  =  je2  dt  +  e 


(A-l) 


where  |f|  <  5  x  10~3  and 
3 


y ■ 


2  2 


(A-2) 


b  a 
\l/3 _ a  ...  —  r 


with  wi  =  (6x)  and  W2  =  («(1  -  x))  .  In  turn,  y  can  be  written  as 
(l  -  x)1/3  (9a  -  l)b213  -  x1/3 (9b  -  l)a2/3 


y  = 


3  a1/2b1/2  Vx2/V/3+(1-x)2/361/3 


(A-3) 
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